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Statistical ensemble formalism of Kim, Mandel and Wolf (J. Opt. Soc. Am. A 
4, 433 (1987)) offers a realistic model for characterizing the effect of stochastic 
non-image forming optical media on the state of polarization of transmitted 
light. With suitable choice of the Jones ensemble, various Mueller transforma- 
tions - some of which have been unknown so far - are deduced. It is observed 
that the ensemble approach is formally identical to the positive operator val- 
ued measures (POVM) on the quantum density matrix. This observation, in 
combination with the recent suggestion by Ahnert and Payne (Phys. Rev. A 
71, 012330, (2005)) - in the context of generalized quantum measurement on 
single photon polarization states - that linear optics elements can be employed 
in setting up all possible POVMs, enables us to propose a way of realizing 
different types of Mueller devices. 
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1. Introduction 

The intensity and polarization of a beam of light passing through an isolated optical de- 
vice undergoes a linear transformation. But this is an ideal situation because, in general, 
the optical system is embedded in some media such as atmosphere or other ambient mate- 
rial, which further modifies the polarization properties of the light beam passing through 
it. A statistical ensemble model describing random linear optical media was formulated 
two decades ago by Kim, Mandel and Wolf [1], but is not examined in any detail in the 
literature, to the best of our knowledge. The purpose of the present paper is to pursue 
this avenue in a new way arising from the realization of a relationship, presented here, with 
the positive operator valued measures (POVM) of quantum measurement theory. This is 
because the transformation of the polarization states of a light beam propagating through 
an ensemble of deterministic optical devices exhibits a structural similarity with the POVM 
transformation of quantum density matrices. This connection motivates, in view of the re- 
cent interest in the implementations of POVMs on single photon density matrix employing 
linear optics elements [2], identification of experimental schemes to realize various kinds of 
Muller transformations. The properties of the transformation of the polarization states of 
light form a much studied topic in literature [3 - 17]. Thus the power of the ensemble 
approach becomes evident in elucidating the known optical devices as well as some hitherto 
unknown types [17], which had remained only a mathematical possibility. 

The contents of this paper are organized as follows. In Sec. 2, a concise formulation of 
the Jones and Mueller matrix theory, along with a summary of main results of Gopala Rao 
et al. [17] is given. Based on the approach of Kim, Mandel and Wolf [1] suitable Jones en- 
sembles, corresponding to various types of Mueller transformations are identified in Sec. 3. 
In Sec. 4, a structural equivalence between Jones ensemble and POVMs of quantum mea- 
surement theory is established. Following the linear optics scheme of Ahnert and Payne [2] 
for the implementation of POVMs on single photon density matrix, experimental setup for 
reahzing Mueller matrices of types I and II are suggested in Sec. 5. The final section has 
some concluding remarks. 
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2. Brief summary of known results on the Jones and the Mueller formalism. 



Following the standard procedure, let Ei and E2, defined here as a column matrix E = 

(eA 

, denote two components of the transverse electric field vector associated with a light 

beam. The coherency matrix (or the polarization matrix) of the light beam is a positive 
semidefinite 2x2 hermitian matrix defined by, 



C = (E®E^). 

Expressing this in terms of the standard Pauli matrices cri = 

, we have 



(1) 



1 

1 
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and the unit matrix (Tq = 


.0 



1 



^ _ 1 ^ ^ _ 1 / So + S3 Si - iS2 ^ 
^ i=0 ^ y Si + tS2 So - S3 J 

The physical significance of the quantities arising here are 



50 = Tr (Cctq) = Intensity of the beam 

51 = Tr (C(7i) = Components of Polarization vector s of the beam 



(2) 



(3) 



Thus the coherency matrix completely specifies the physical properties of the light beam. 

The four-vector S = defined by Eq. (3) is the well known Stokes vector, which 

\^ J 

represents the state of polarization of the light beam. Because C is hermitian, the Stokes 
vector is real. The positive semidefiniteness of C implies that the Stokes vector must satisfy 
the properties 

So >0, si - |s|^ > (4) 

A 2x2 complex matrix J, called the Jones matrix, represents the so-called deterministic 
optical device [18] or medium. When a light beam represented by E passes through such a 
medium, the transformed light beam is given by E' = JE. Correspondingly, the coherency 
matrix C transforms as 

C = JCJ^ (5) 

(Here is the hermitian conjugate of J.) 
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Alternatively, instead of the 2x2 matrix transformation of the coherency matrix, as 
given by Eq. (5), a transformation 

S' = MS (6) 

of the four componental Stokes column S through a real 4x4 matrix M, called the Mueller 
matrix, is found be more useful [18]. 
Using Eq. (3) and Eq. (5 we have, 

4 = Tr(CVO = Tr(JCJt a^) = ^ ^Tr(JtaiJ(7,•)s,■ 
which leads to the well-known relationship [1] 

between the elements of a Jones matrix and that of corresponding Mueller matrix. 

But in the case where medium cannot be represented by a Jones matrix, it is not possible 
to characterize the change in the state of polarization of the hght beam through Eq. (5). 
In such a situation, Mueller formalism provides a general approach for the polarization 
transformation of the light beam. The Mueller matrix M is said to be non-deterministic 
when it has no corresponding Jones characterization. 

Mathematically, a Mueller device can be represented by any 4x4 matrix such that the 
Stokes parameters of the outgoing hght beam satisfy the physical constraint Eq. (4) . In other 
words, a Mueller matrix is any 4x4 real matrix that transforms a Stokes vector into another 
Stokes vector. There are many aspects of the relationships between these two formulations 
of the polarization optics and a complete characterization of Mueller matrices has been the 
subject matter of Ref. [1, 3-17]. It was Gopala Rao et al. [17] who presented a complete set 
of necessary and sufficient conditions for any 4x4 real matrix to be a Mueller matrix. In so 
doing, they found that there are two algebraic types of Mueller matrices called type I and 
type II; and it has been shown [17] that only a subset of the type-I Mueller matrices - called 
deterministic or pure Mueller matrices - have corresponding Jones characterization. All the 
known polarizing optical devices such as retarders, polarizers, analyzers, optical rotators are 
pure Mueller type and arc well understood. Mueller matrices of the Type II variety arc yet 
to be physically realized and have remained as mere mathematical possibility. For the sake 
of completeness, we present here the characterization as well as categorization of these two 
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types of Mueller matrices as is given in Ref. [17]. This will enable us to show that both 
Type I and II Mueller devices are realizable in an unified manner in terms of the proposed 
ensemble approach [1]. 

I. A 4 X 4 real matrix M is called a type-I Mueller matrix iff 

(i) Moo > 

(n) The G-eigenvalues po, pi, p2, ps of the matrix N=MGM are all real. (Here, M 
stands for the transpose of M; G-eigenvalues are the eigenvalues of the matrix GN, 
with G = diag(l, -1, -1, -1)). 

iii) The largest G-eigenvalue po possesses a time-hke G-eigenvector and the G-eigenspace 
of N contains one time-hke and three space-hke G-eigenvectors. 

II. A 4 X 4 real matrix M is called a type-II Mueller matrix iff 

(i) Moo > 0. 

(n) The G-eigenvalues po, pi, P2, Ps of N=MGM are all real. 

(iii) The largest G-eigenvalue po possesses a null G-eigenvector and the G-eigenspace of N 
contains one null and two space-like G-eigenvectors. 

(iv) If Xo = Bo -|- 61 is the null G-eigenvector of N such that eo is a time-like vector 
with positive zeroth component, ei is a space-like vector G-orthogonal to eo then 
eoNeo > 0. 

Despite the knowledge of these new category of Mueller matrices [15, 17], not much 
attention is paid for realizing the corresponding devices. An experimental arrangement 
involving a parallel combination of deterministic (pure Mueller) optical devices is proposed 
in Ref. [17] for realizing type-II Mueller devices. The physical situations, where the beam of 
light is subjected to the influence of a medium such as atmosphere was addressed in Ref. [1]. 
In the next section, we discuss this ensemble approach for random optical media, proposed 
by Kim, Mandel and Wolf [1] . 
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3. Mueller matrices as ensemble of Jones devices 

Kim et. al. [1] associate a set of probabilities {pe, J2Pe = 1} to describe the stochastic 
medium. Then a Jones device Je associated with each element e of the ensemble gives a 
corresponding coherency matrix C'e = JgCJ^. The ensemble averaged coherency matrix 

Ca,=J2PeC'e = T.PeiJeCJl) (7) 

e e 

then describes the effects of the medium on the beam of light. In a similar fashion, the 
corresponding ensemble of Mueller matrices {Mg} associated with the ensemble of Jones 
matrices {Je} is constructed and its ensemble averaged Mueller matrix is similarly formed 
as M.av — Z^ePeMg. Since a linear combination of Mueller matrices with non-negative 
coefficients is also a Mueller matrix, the ensemble averaged Mueller matrix M^^ is a Mueller 
matrix ^. 

We now turn to the question of constructing an appropriate ensemble designed to describe 
a given physical situation. The simplest example of an ensemble is one where the elements 
are chosen entirely randomly, i.e., the system is described by a chaotic ensemble where the 
probabilities are all equal, Pe = ^, where n denotes the number of elements in the ensemble. 
The coherency matrix Co„ of the light beam passing through such a chaotic assembly is just 
an arithmetic average of the coherency matrices C'e = JeC J^ and hence 

1 " 

C„.= -EJeCJi (8) 

More general models can be constructed depending on the medium for the propagation of 
the beam of light. For example, one may employ various types of filters or solid state systems 
through which the light passes; the assignment of the Jones matrices and the corresponding 
probabilities will then differ depending on the weights placed on these elements. 

Restricting ourselves to an ensemble consisting of only two Jones devices which occur 
with equal probability pi = 1/2, p2 = 1/2, we have found out that the resultant Mueller 
matrices can either be deterministic or non-deterministic. We give in the foregoing (see Table 
I) some examples of Mueller matrices corresponding to different choices of Jones matrices 
in an ensemble Jg, e = 1, 2, for some representative cases. This will also serve to show the 

^ This is because, each Mueller matrix Me transforms an initial Stokes vector into a final Stokes vector and 
a linear combination of Stokes vectors with non-negative coefficients Pe is again a Stokes vector. 
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generality of the ensemble procedure in capturing the physical realizations for the Mueller 
devices discussed in Ref. [17]. 



Table 1. Abieller inatrices resulting from 2-eleiiierit Jones ensemble. 



J2 



M = piMi +P2M2, 

Pl=P2 = h- 



Type of M 



'■is 



1 1 

l + i - 



1 

V6 



1 1 ■ 



/ 3 \ 
0-122 

2-12 

\ 2 2 -1 / 



Pure Mueller 





/ 1 \ 

10 



\ / 



Type-I 



3. ^ 



1 

1 



1 

■/2 



1 
-1 



/ 1 \ 





\ -1 / 



Type-I 



4- 4- 

\/3 



1 

1 




/ 3 \ 

0-102 
0-10 
V 2 -1 / 



Type-I 



1 1 -i 
l + i -1 



1 



1 -i 
i -1 



10 



/ 5 \ 

0-124 

2-32 

\ 4 2 -1 / 



Type-I 





/ 1 -1 \ 



\ 000/ 



Type-II 



1 1 
1 -1 



/ 2 1 \ 

10 

10 

V 1 / 



Type-II 



In Table I, the Jones matrices chosen are so as to give pure Mueller (deterministic) 
and non-deterministic type-I, type-II matrices respectively. We observe that an assembly 
of Jones matrices can result in a pure Mueller matrix if and only if all elements of the 
assembly correspond to the same optical device. This is because, with all Je's are same, 
a transformation of the form Cav = J2ePe{JeC3l) is equivalent to a transformation of the 
Stokes vector S through a Mueller matrix M^^ = Y^p^Me — Mpure- When the medium is 
represented by a pure Mueller matrix, the outgoing light beam will have the same degree of 
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polarization as the incoming light beam. In fact, pure Mueller matrix is the simplest among 
type-I Mueller matrices. Not all type-I Mueller matrices preserve the degree of polarization 
of the incident light beam. To see this, note that the type-I matrix of example 2 (see Table 
I) converts any incident light beam into a linearly polarized light beam; the other three 
type-I matrices (examples 3 to 5) transform completely polarized light beams into partially 
polarized light beams. Similarly, type-II Mueller matrices do not, in general, preserve the 
degree of polarization of the incident hght beam. It may be seen that the type-II Mueller 
matrix of example 7 is a depolarizer matrix, since it converts any incident light beam into 
an unpolarized light beam. 

Though one cannot a priori state which choices of Jones matrices result in type-I or type- 
II, it is interesting to observe that all types of Mueller matrices result - even in 2-element 
ensembles. It is not difficult to conclude that an ensemble, with more Jones devices and 
with different weight factors, can give rise to a variety of Mueller matrices of all possible 
algebraic types. It would certainly be interesting to physically realize such systems. 

In the following section, a connection between the ensemble approach for optical devices 
and the POVMs of quantum measurement theory is established. 

4. A connection to Positive Operator Valued Measures 

We will now show that the phenomenology of the ensemble construction of Kim, Mandel 
and Wolf [1] described above has a fundamental theoretical underpinning, if we make a 
formal identification of the coherency matrix with the density matrix description of the 
subsystem of a composite quantum system. The coherency matrix defined by Eqs. (1) 
and (2) resembles a quantum density matrix in that both describe a physical system by 
a hermitian, trace-class, and positive semi-definite matrix. While the quantum density 
matrix has unit trace, the coherency matrix has intensity of the beam as the value of the 
trace. The Jones matrix transformation is a general transformation of the coherency matrix, 
which preserves its hermiticity and positive semi-definiteness - but changes the values of the 
elements of the coherency matrix. The most general transformation of the density matrix 
p, which preserves its hermiticity, positive semi-definiteness and also the unit trace is the 
positive operator valued measures (POVM) [19]: 

n n 

p' = EV.pVl; EVIV, = I (9) 

i=i i=i 
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where T^'s are general matrices and I is the unit element in the Hilbert space. More generally, 
one could relax the condition of preservation of the unit trace of the density matrix by 
examining the possibility of a contracting transformation, where the unit matrix condition 
on the POVM operators is replaced by an inequality. 

This mathematical theorem has a physical basis in the Kraus operator formalism [19] 
when we consider the Hamiltonian description of a composite interacting system A, B 
described by a density matrix B) and deduce the subsystem density matrix of A given 
by, p{A) — T^B p{A, B). In this case, the Kraus operators are the explicit expressions of 
the POVM operators and contain the effects of interaction between the systems A and B 
in the description of the subsystem A. It is thus clear that the phenomenology of Rcf. [1] 
has a correspondence with the Kraus formulation and the POVM theory. In order to make 
this association complete, we compare Eq. (9) with the expression given by Eq. (7). Apart 
from a phase factor, the Kraus operators {Vj}, associated with POVMs, may be related to 
the Jones assembly {Jj}, chosen in the form 

V, = J„ vjVi = Jl Ji (10) 

i=l i=l 

In the construction of the Table I presented earlier, a simple model was proposed where all 
probabilities were chosen to be equal and the condition on the sum over the Jones matrix 
combinations was set equal to unit matrix. In such cases, the intensity of the beam gets 
reduced by 1/n and the polarization properties of the beam gets changed as was described 
earlier. With this identification, we have provided here an important interpretation and 
meaning to the phenomenology of the ensemble approach of Kim et al. [1] . 

Recently Ahnert and Payne [2] proposed an experimental scheme to implement all possible 
POVMs on single photon polarization states using linear optical elements. In view of the 
connection between the ensemble formalism for Jones and Mueller matrices with the POVMs, 
a possible experimental realization of the two types of Mueller matrices is suggested in the 
next section. 

5. Possible experimental realization of types I and II Mueller matrices. 

We first observe that the density matrix of a single photon polarization state, 

p = phh\H) {H\ + pHv\H) {V\ + p*Hv\V) {H\ + pvv\V) {V\ (11) 
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is nothing but the coherency matrix of the photon [20] 

C = ( ^^^^"^ ^^^^^^ I , (12) 

where a^^ and ky are the creation operators of the polarization states of the single photon; 
{\H), \V)} denote the transverse orthogonal polarization states of photon. This is seen 
explicitly by noting that the average values of the Stokes operators are obtained as, 

So = (So) = ((aj^a// + a|.ay)) = phh + Pvv = Tr(p), 

si = (Si) = ((ajjay + alkn)) = Phv + P*hv ^ Tr(p(Ji), 

52 = (S2) = i {{a^eLH - aj^ay)) = i {pnv - Phv) = Tr(pcr2), 

53 = (S3) = ((a]^a// - a\.av)) = Phh - Pvv = Tr(p(73). (13) 

Hence the proposed setup [2], involving only linear optics elements such as polarizing beam 
splitters, rotators and phase shifters, that promises to implement all possible POVMs on a 
single photon polarization state leads to all possible ensemble realizations for the Mueller 
matrices. More specifically, this provides a general experimental scheme to realize varieties 
of Mueller matrices - including the hitherto unreported type-II Mueller matrices. We briefly 
describe the scheme proposed in Ref. [2] and illustrate, by way of examples, how it leads to 
both type-I and type-II Mueller matrices. 

In Ref. [2], a module corresponds to an arrangement having polarization beam splitters, 
polarization rotators, phase shifters and unitary operators. For an n element POVM, a 
setup involving n — 1 modules are needed. That means, a single module is enough for a 2 
element POVM; a setup involving two modules is required for a 3 element POVM and so on. 
We describe two, three element POVMs by specifying the optical elements in the respective 
modules and by specifying the corresponding Kraus operators in terms of these elements. 

For any two operator POVM, the Kraus operators Vi, V2 are given by Vi = U'DiU 
and V2 = U"D2U. Here U, U', U" are the three unitary operators in a single module. 
Denoting 6, as the angles of rotation of the two variable polarization rotators and 7, ^, 
the angles of the two variable phase shifters in the module, the diagonal matrices Di, D2 
are given by. 



e^Tcos^ \ I e'^sind ^ 

Di = I ' = 

COS0 / \ sin0 



(14) 
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The POVM elements 



Fi = VjVi = U^dIDiU, F2 = V^V2 = U^D^DsU 



(15) 



satisfy the condition J2i=i,2 Fj = Fi + F2 = I. 

For any three operator POVM, the Kraus operators are given by 

Vi = UiDiUi, 

V2 = UiiDnU,,Ui'DiUi, 

V3 = Ui'iD'„U„U'/DiUi, 



(16) 



Here, the diagonal D matrices are 
Di = 



e'Ticos^i 
cos (f)i 



and 



D 



II 



e*T"cos^ii 
cos 011 



e*^i sin 61 
sin 01 



^ e*^" sin 611 



(17) 



II 



V 








(18) 



(^i, 0i), (71, ^i) are respectively the pair of angles corresponding to variable polarization 
rotators and variable phase shifters in the first module. Similarly, (^n, <Pii), {in, Cii) the 
pairs of angles corresponding to variable polarization rotators and variable phase shifters 
respectively in the second module. Ui, Uj, U'/ arc the unitary operators used in the first 
module and Un, Ujj, Ujj are the unitary operators used in the second module. (Notice 
that all the unitary operators in the above schemes are arbitrary and a particular choice of 
the associated unitary operators gives rise to a different experimental arrangement). The 
extension of this scheme to n operator POVM involving n-1 modules is quite similar and is 
given in [2]. 

We had identified, in Sec. 3, that an ensemble average of Jones devices will lead to all 
possible types of Mueller matrices, some examples of which are given in Table 1. We now 
show that the experimental set up proposed in Ref. [2] can also be used to realize varieties 
of Mueller devices. To substantiate our claim, we identify here the hnear optical elements 
needed in the single module set up of Ahnert and Payne [2], which lead to the physical 
reahzation of two typical Mueller matrices given in Table 1. 
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To obtain the type-I Mueller matrix M= 



1 











-1 



of example 3 (see Table I) , we 



use U = I, U' 







and U" = 


: 











1 
-1 



as the required unitary Jones devices and 



both the variable polarization rotators are set with their rotation angles 9—(f) — n/A. There 



is no need of phase shifter devices in this case i.e, 7=^=0. 

/ 



Similarly for the type-II Mueller matrix M= 



2 1 
10 
10 
10 



of example 7, we find that 



1 1 
1 -1 



and U" = ^ 



1 i 
i 1 



are the required unitary Jones devices. The 



rotation angles of the variable polarization rotators are, as in the earlier case, 9—(f) — n/A 
and there is no need of phase shifter devices i.e., 7=^=0. Notice that in both the above 
examples the unitary operators U', U" correspond to linear and circular retarders [18]. 

These two examples illustrate that the experimental set up given in Rcf. [2] may be 
utilized to realize the required non-determinisitc Mueller devices. In fact Mueller matri- 
ces corresponding to an ensemble with more than two Jones devices may also be realized 
by employing larger number of modules as given in the experimental scheme proposed by 
Ref. [2]. 

6. Conclusion 

We have established here a connection between the phenomenological ensemble approach [1] 
for the coherency matrix and the POVM transformation of quantum density matrix. This 
opens up a fresh avenue to physically realize types I and II of the Mueller matrix classification 
of Ref. [17]. We have also given experimental setup to implement Mueller transformations 
corresponding to ensemble average of Jones devices by employing the POVM scheme on the 
single photon density matrix suggested in Ref. [2], in the context of quantum measurement 
theory. It is gratifying to note that two decades after the introduction of the ensemble 
approach, which had remained obscure and only received passing reference in textbooks such 
as [20] , its value is revealed in this paper through its connection with the new developments 
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in quantum measurement theory. We plan on exploring further the POVM transformation 
in the description of quantum polarization optics. 
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